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ABSTRACT 

The thermospheric wind circulation, excited during magnetic storms 
presujnably by joule heating within the auroral zone, is shown to be 
an effective mechanism for removing atomic oxygen at high lati- 
tudes. Wind induced variations in O exceed the temperature effects 
up to 250 km. The calculated deplation is most pronounced at around 
180 km where the density can decrease by as much as a factor of 
two, consistent with the observed storm time variations in the 
ionosphere. At higher altitudes, this effect is canceled by the ther- 
mal expansion in atomic oxygen thus explaining the negligible re- 
sponse in the concentration of this atmospheric constituent under 
distributed conditions when increased by as much as a factor of 
ten. 
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MAGNETIC STORM EFFECTS IN THE NEUTRAL COMPOSmCN 


by 

H. G. Mayr and H. Volland 


I. INTRODUCTION 

Reber et al. 1970 and Taeusch and Carignan 1970, have presented mass spec- 
trometer measurements of N 2 and O from OGO 6 during magnetic storms. These 
data show at heights of typically 400 km the following characteristics: 

1. A gradual increase in the N^ concentration from low to high latitudes 
(up to a factor of 10 for ^ 200), 

2. An increase in the exospheric temperature (from 1000° up to 1400°) at 
high altitudes which was infered from the variations in N 

3. Negligible variations in O with some tendency to decrease at high latitudes 
where Nj showed significant enhancements. 

The obvious implication is that during magnetic storms a considerable 
amount of energy is deposited in the auroral zone presumably below 200 km 
thus causing a thermal expansion in Nj which is the major constituent in the 
lower thermosphere. 

Associated with this temperature enhancement, O should increase similarly 
if, as present models assume, the composition did not vary at the turbopause 
level. This is contrary to the observed behavior in atomic oxygen which, if 
anything, tends to decrease in the region of enhanced temperatures. 


Thus we suggest that the magnetic storm heating produces changes in the 
neutral composition and in particular a decrease of the relative concentration of 
O withinthe lower thermosphere. In this paper, we shall discuss a model that 
can explain such variations in terms of diffusion processes associated with the 
thermospheric circulation that is excited in an altitude range where joule dissipa- 
tion (Cole, 1970) is most likely the predominant heat source. 

II. THEORY 

The thermospheric winds up to 250 km, where Nj is the major constituent, 
exert drag forces on the minor species Oin this region. The diffusion velocites 
induced by this interaction upset the static equilibrium for atomic oxygen thus 
producing variations in the relative concentration of O to N 2 . This process has 
been studied for seasonal variations in the annual component (Kellogg, 1961) 
and semiannual component (Mayr and Volland, 1970). 

Magnetic storms have time scales of typically one day which are quite 
different from the periods of half a year in the seasonal variations, and this 
reflects significantly upon the atmospheric properties. One can estimate (e.g. 
at 120 km) the characteristic times for diffusion of O through 

, , 

-x. 6 days 

D D 

and for the mass transport due to horizontal winds 



where H(10 km) is the density scale height, D(6 x 10^) the diffusion coefficient, 
V (10 m/sec) the wind velocity and R the earth’s radius. 
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Seasonal variations in the wind field and in the diffusion transport can easily 
be excited even at the turbopause level. The response of the neutral atmosphere 
is almost instantaneous on this time scale thus the continuity equations can be 
treated as time independent. As a consequence, a unique relation evolves 
between the vertical and horizontal wind velocities with upward winds reducing 
atomic oxygen and downward winds enhancing it (Mayr and Volland, 1970). 

In contrast, magnetic storm variations with a time scale of about one day, 
comparable to the response time of the atmosphere, cannot be excited readily at 
all altitudes. The response time of the wind field and the phase difference be- 
tween vertical and horizontal winds increase significantly at low altitudes. Tho 
relatively long diffusion time causes further delay in the response of the com- 
position to the vertical and horizontal wind fields; so it is obvious that the tran- 
sient aspects of the magnetic storm variations must be important. 

We shall therefore address ourselves to the following questions: can the 
wind fields during magnetic storms produce variations in the composition that 
are consistent with the observations and b) what inferences can be made on the 
energy source from the time response of the composition? 


The Model of Volland and Mayr (1971) 

The bulk motion of up to 250 km where this species dominates constitute 
the wind field that affects theO concentration. As a major constituent, Nj is 
not significantly influenced by diffusion. Furthermore, the energy sources that 
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determine essentially these dynamics are within the Nj region up to 250 km. For 
these reasons, the wind field can be computed without considering diffusion pro- 
cesses. However, the wind fields at higher altitudes where O is the major con- 
stituent is of course affected by these processes. This aspect will be discussed 
with some of its consequences in a subsequent paper. 

Volland and Mayr (1971 a,b) have treated the dynamic of the neutral atmosphere 
in response to auroral heating. Although the calculations cover an altitude range 
up to 400 km diffusion has not been considered. For reasons outlined before, 
however, their results are completely valid within the region of the thermo- 
sphere and thus can be adopted to study the diffusion properties of 0. A brief 
summary of this model is therefore here appropriate. 

Assuming an energy source that is confined to the auroral zone, Volland 
and Mayr expanded this input into a series of spherical harmonics 

(Pq + 3.66 f 2.22 - 2.65 Pg ~ 6,91. Pg • • . . ) (1) 

As expected, this series converges slowly. The model shows however, that 
the response of the neutral atmosphere to the various energy components is such 
that the higher wave numbers are rapidly damped. This is illustrated in Fig. 1 
(corresponding to Fig. 2 in their paper) which shows the transfer function 
that relates density and energy (Q^) variations over a wide frequency 
range through 

P - G Q , (2) 

thus describing the efficiency for the excitation of the various components. As 
is evident from this figure, the efficiency decreases with the power of l/n^ for 
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frequencies that correspond to magnetic storm periods of typically one day 
(o! = 7.3 X 10"** /sec). This clearly implies that only the first two harmonics 
Pq and in the energy expansion are most effective in exciting density variations, 
a fact that is supported by measurements that show a strong worldwide com- 
ponent (Pq) in the magnetic storm response of the thermosphere and a gradual 
increase toward higher latitudes (Pj ) in the Nj concentration (Taeusch and 
Carignan, 1971). Thus we shall confine ourselves in this paper to treat only the 




whero 


Iv Boltzimoim constant 


q ~ photo production rate due to dissociation of O 
= rate of three body recombination (0 + 0 + Nj Oj + Nj) 

= rate of three body recombination (O + Oj + Nj O3 + N3) 

Oh OjS, IN2 ' ~ number densities for O, O2 and Nj 
Vqt* Vqp, = transport velocities for O in the radial (r), colatitudinal 
(i‘) and longitudinal (A.) direction. 

Vf. V^, ~ wind components for the three directions 

D = molecular diffusion coefficient for 0 through Oj and 
K = Eddy diffusion coefficient 
t = universal time. 

m,m = mass of atomic oxygen^ mean molecular mass 
g » grxivitatlonal acceleration 

In the latitudinal and longitudinal directions the drag terms tO] ( - v^) 

and iOj (Vqx ■ Vx ) dominate the lateral momentum transfer thus reducing to 
the simplified equations ( 4 ) and ( 5 ). The inertia terms B/ 3 t (Vq^, have 

been neglected in Equations 3 through 5 since they are negligible when compared 
with accelerations due to the pressure gradients. 

In order to solve these equations, a number of simplifying assumptions were made: 

1 . The local time and longitude dependences are neglected 

2. The variables are separated into storm time t, colatitude and altitude 
r‘ components in the form 

[0] = [0]o(r)+ [OljCr) (6) 

for atomic oxygen, and 
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F" 101 (V„, -V/) - 

•• F(,(i-) 1 Fj(r) PjC'O (7) 

for the diffusion flux. The subscripts 2 refer to the second spherical harmonic 

PjCfO ^ (3 cos^ 0 - 1). 

We choose the predominant Fourier component for the magnetic storm variations 
thus adopt a frequency a with a period of typically one day. 

3. From the model of Volland and Mayr (1971, b), which applies full wave 

calculations to the atmospheric dynamics, the appropriate forms for the 
wind and temperature fields are adopted. 

( 10 ) 

Fg =Ffl,(r) sin^)P|(e) e’*“' (11) 

with 

Pj((9) - cos{B), 

and 

T = To(r) ^ TjCr) Pj (12) 

4. Assuming that higher order terms in the frequency and latitude expansions 
can be neglected, perturbation theory is applied which leads in a straight forward 
manner to a set of first order differential equations 

~ (Fq) ^ aj [0]2 [Nj] ^ aj [0]^ [Oj] [Nj] - 2q [Oj] = 0 (13) 
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where N is the atmospheric density. 
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The equations (13) and (14) describe the steady state distribution for atomic 
0 x 3 gen without being affected by winds and temperature variations. As these 
equations do not depend on the other variables, they can be solved independently. 
Such a solution has been described e.g. in Mayr and Volland (1970) and it is here 
adopted as input. 

The equations (15) and (16), which are complex and represent therefore 
four equations, are then solved by means of a trial and error scheme with the 
following boundary conditions: 

a. At heights above 250 km,0 becomes the major constituent, therefore 
diffusion is negligible 


f = 0 at 250 km 

b. The diffusion time increases at altitudes below 100 km to four days, thus 
the temporal variations must be small there. To account for that, it is assumed 
arbitrarily that 

P = 0 at 90 km. 

This implies that the solution is meaningless at this altitude level. However , at 
heights above 100 km, the solution is very insensitive to this artificial constraint 
thus the solution can be considered as unique there. 

The er lations (15) and (16) have the important characteristic that the complex 
variables f and p depend linearily on the temperature and wind fields V ^2 
and V^j . This implies that the effects induced by temperature and wind varia- 
tions can be separated thus providing a means of identifying their relative 
significance. 
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Considering furthermore that the temperature and wind fields from the model 
of Volland and Mayr (1971 b) also depend linearily on the energy input, a linear 
relationship between variations in the composition and heat source evolves 
which greatly simplifies the discussion. 

III. DISCUSSION 

Two significantly different heating mechanisms have been proposed for 
geomagnetic storms in the Auroral zone: 

a. Heating due to particle precipitation, a process in which recombination, 
electron heating and vibrational excitation of N 2 produce a heat input that peaks 

at about 100 km or above, depending upon the energy of the precipitating electrons. 
(Rees, private communication.) 

b. Joule heating due to electric fields that are generated in the magnetosphere 
(Cole, 1970). Particle precipitation contributes to this heat source by enhancing 
the conductivity through impact ionization. This energy input peaks in the E 
region presumably at about 150 km (Cole, 1970). 

It seems obvious that both mechanisms should be active. The difference in 
the height distribution for both heat inputs, however, suggests significantly dif- 
ferent responses in the thermosphere dynamics. Thus the question arises whether 
the observed effects in the composition can be tied (through our models) to one 
particular heat source ? 

Two energy models, shown in Fig. 2 are adopted. The form of the distribu- 
tions is the same , however the energy peaks are placed at 100 and 150 km. As 
we deal here with a linear theory, the absolute magnitudes of the heat input is 
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insignificant at this point. The relative magnitudes of the and components 
are chosen in accordance with the spherical expansion (1) 

Q 2 ^ 3.66 Qq. 

With these inputs, the temperature and wind fields are computed from the three 
dimensional model of Volland and Mayr (1971, b) for the predominant frequency 
of a typical magnetic storm {co = 2tt/ 1 day). As we restrict ourselves here to the 
dominating frequency we make a very rough approximation to the impulse type 
heat input in the Auroral zone. For this reason our present analysis cannot ac- 
count for details in the time evolution of the storm time dynamics. 

The results show that the uniform (one dimensional) expansion for the Qq 
component produces vertical (only) wind velocities that are one order of magni- 
tude smaller than the vertical wind velocities for the non-uniform and two 
dimensional Pj component, while the temperatuer amplitude is somewhat larger 
for the Qq term. As we are primarily interested in wind- diffusion effects com- 
peting with the temperature induced variations in the composition, we conclude 
that the Q q component-equivalent to a one-dimensional approximation - does not 
contribute appreciably to the dynamics and diffusion in the thermosphere. 
Therefore, we can disregard the term. 

Wind Effects 

Fig. 3 shows the magnitudes of the vertical and meridional wind velocities 
and temperature fields computed from the model of Volland and Mayr (1971, b) 
for the component. Figs 3a and 3b represent the solutions for energy inputs 
that peak at 100 km and 150 km respectively (see Fig. 2). 
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The wind direction is upward (positive) for the vertical wind field in Figs. 

3a and 3b at the time of the maximum Auroral heat input. At the same time the 
meridional wind direction is equatorward (positive) at heights of the energy maxi- 
mum and above (that is above 100 km for Fig. 3a and above 130 km for Fig. 3b). 
Below 130 km the winds blow toward the poles for case b. 

With these inputs, the diffusion equations (15) and (16) are solved. Fig. 4a 
describes amplitudes and phase of the variations in atomic oxygen for the vertical 
and meridional wind components and for the combination of both. These varia- 
tions are excited by a heat input that beaks at 100 Ion. Fig. 4b describes these 

I 

parameters for the energy maximum at 150 km. The time of the density maxi- 
ma is plotted with respect to the time of the maximum auroral heat input. 

The most interesting feature in the comparison between both solutions — that 
show the combined wind effect (solid lines in Fig. 4) for these energy 
models — is the phase relation. In both cases, the phase is very similar below 
120 km. Above 140 km, however, a phase difference of nearly 12 hours deve lopes 
thus indicating that the effects go in opposite directions. As is evident, the 
maximum in atomic oxygen occurs essentially at the time of the maximum auroral 
heat input (Fig. 4a) in the case where the wind field is excited at 100 km. For the 
150 km heat input, this maximum occurs 12 hours (180°) earlier nd this implies 
that a minimum developes in atomic oxygen during the time of magnetic storms. 

The reasons for this difference in the response of the neutral composition are 
somewhat complex. The vertical and meridional wind velocities enter the diffu- 
sion equation (15) in the form 
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> 0 (sou r ce) 


V 


r2 




<' 0 (si nk) 


(20) 


where 


H, 



(oxygen scale height) 


and 



r2 



(scale height of the vertical winds) 


In Fig. 3 the meridional velocity is positive (except below 130 
km in Fig. 3b and there the wind velocity decreases rapidly). For this reason, 
this wind component acts as a sink which removes atomic oxygen at high altitudes. 
This is reflected both in Figs. 4a and 4b where the component (AO)v is shifted 
to earlier hours thus producing a minimum in atomic oxygen at or after the 
time of the maximum energy input. At heights below 140 km, the phase is shifted 
progressively to later hours at lower heights thus reflecting the increasing re- 
sponse time in the lower thermosphere. 

The effects of the vertical wind velocity depend on its scale height, 

(see Eq. 20). When the vertical wind velocity increases more rapidly than 
atomic oxygen is decreasing with height the wind effectively removes oxygen; 
otherwise it supplies atomic oxygen. In this respect, particularly the wind fields - 
excited at 100 and 150 km - are very different. 

In the regions around the energy source and up to about 130 km atomic oxygen 
is removed to be supplied over the entire altitude range between 130 and 250 km 
(for Fig. 4a). Above 130 km the vertical wind constitutes a source that accumulates 
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oxygen during storm times. This is reflected in both, the relatively large amplitude 
of (/\0) (which surpasses the component (AO)„ ) and in the occurance of the 

density maximum after the time of the maximum energy input. Both wind com- 
ponents (AO)v and (AO)y are out of phase by almost 180° thus opposing each 
other. The predominance of the (AO)y ^ term therefore produces a slight en- 
hancement in atomic oxygen during storm time (see the phase and amplitude of 
(A0)^.|^^ in Fig. 4a). This is opposite to the variations in the composition that 
are inferred from storm time measurements. Below 130 km the steep gradient 
in causes a depletion of atomic oxygen which is reflected in the decrease 
of the amplitude (AO)„ and in the shift of its phase to earlier hours. The 

'^r2 

superposition of (AO)y and (AO)^ produces then the structure in (AO)^^^^j 
within the lower thermosphere. 

In the case where the energy source peaks at 150 km, the vertical wind 
component removes atomic oxygen up to 170 km and supplies it over an altitude 
range of only 50 km above that height. The net effect is a relatively small ac- 
cumulation of O above 150 km which is reflected in amplitude and phase of (AO)^ 

1-2 

(illustrated in Fig. 4b). In contrast to Fig. 4a, the effect of the vertical wind 
field is here considerable smaller than that of the meridional wind velocity. 
Consequently, the depletion ofO due to meridional winds dominates (see (AO)^-^^ 
in Fig. 4b) and O decreases during storm time. 

Below 170 km, the effects are similar to the ones below 130 km which were 
described in Fig. 4a. 

Finally, a very important feature in the comparison between both models 
is related to the values of the phase. For the energy source at 100 km, the 
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vertical wind velocity produces a peak and the meridional wind a minimum in O 
five hours after the storm above 140 Ion. With the energy input at 150 km, the 
phase difference is only an hour, thus contributing to the depletion ofO at the time 
of the magnetic storm. This difference simply reflects the decrease of the response 
time in the neutral atmosphere at higher altitudes. 

Comparison Between Temperature and Wind Effects 

In Figs. 5a and 5b, we show the variations inO due to the thermal expansion 
(AO)^ (computed from Eq. 15 and 16), and compare them with the wind induced 
components that were discussed before. Furthermore, we plotted the variations 
in Nj. In order to cover the altitude range shown in these figures, we have 
assumed that the wind effects are negligible above 250 km. Therefore phase 
and amplitude remain constant. (Ao)^ was deduced above 250 km bj 

considering the diffusion equations that express the state of hydrostatic equilibrium 
there. The distributions in N 2 were adopted from the total atmospheric density up 
to 250 km in the three dimensional model of Volland and Mayr (1971, b) and they 
were extrapolated above 250 km under the assumption of diffusive equilibrium. 

Figs. 5a and 5b show the two models for energy maxima at 100 and 150 km. 

It is apparent from these results that the magnitudes of the wind induced 
variations in atomic oxygen are comparable or even exceed the temperature 
components. 

Above 200 km, the thermal expansion causes in both models a maximum in O 
at the time of the maximum auroral heat input. This is in phase with the time 
variation of N 2 . Below 200 km, the phases of the temperature components shift 
in accordance with the phase and height gradients of the temperature field. 
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The superpositions oi the temperature and wind components in the variations 
oi‘ atomic oxygen are shown as thick lines with the annotation AO. 

For the energy input at 100 Ian, the amplitude increases with height after 
showing a structure below 140 km which is somewhat similar to the varia- 
tions in N 2 .* 

The phase remains almost constant over the entire altitude range showing 
atomic oxygen to peak in phase with at the time of the maximum heat input. 
Accordingly, we conclude that this energy model produces a composition effect 
that is contrary to the observations. 

The results in Fig. 5b, for an energy maximum at 150 km are completely 
different. The amplitude peaks at 180 km with the phase at- 12 hours. Thus 
atomic oxygen decreases and a minimum deve lopes at the time of the storm. As 
is apparent in Fig. 5b, this effect is essentially induced by atmospheric winds. 

At higher altitudes, the temperature induced variations become more and more 
significant. Since temperature and wind effects go in opposite directions (see 
the phase difference of nearly 12 hours = 180°) they cancel each other and 
thus produce the decrease in the amplitude and the gradual shift in the phase 
from -12 to -6 hours between 220 and 340 km. Above that height, the thermal 
expansion dominates the variations in O, thus causing some recovery in the 
amplitude and a further shift in the phase toward the phase of the temperature 
component at the time of the maximum energy input. As a result, the storm 
time increase at higher altitudes (-^400 km) is significantly reduced which is 
in agreement with the observations by Reber, et. al. (1970) and Taeusch and 
Carignan (1970). 

*This structure is of no significance; it reflects only the somewhat artificial boundary conditions 
in the wind model at 100 km which are damped out at higher altitudes. 
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Figure 6 summarizes these results. In altitude intervals of 100 km, we 


show the ratios 


0]^ - 2o0 Pj 
_ 


and 


NjJq + 2AN2P2 


for t = 0 (at the time of the maximum energy input), which reflect the relative 
variations inO and during a magnetic storm. 


It is apparent from this figure that N 2 increases toward higher latitudes 
and at all alti*^udes thus showing the effect of the thermal expansion. In contrast, 

O is shown to decrease by about a factor of two at lower heights and at mid to 
high latitudes. This is consistent with the ionospheric analysis of Chandra and 
Herman (1969) who inferred a corresponding decrease of theO/N 2 ratio during 
storms. At 300 km, no variation is apparent. At 400 km, O increases only slightly 
while N 2 increases by about a factor of 10 from low to high latitudes; both features 
inO and N 2 are in agreement with the observations of Taeusch and Carignan, 1970. 


To reproduce the observations , a peak energy input of 

Q 2 = 0.49 er g/cm^ sec 

was required. Correspondingly, the global energy component is according to (1) 

Qq = 0,13 er g/cm^ sec 

This energy input is comparable with the EUV heat source from the solar radiation. 
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IV, SUMMARY AND CONCLUSION 


Fig. 7 summari?iOS iu schematic iorm our niodcl for the storm time varia- 
tions in the neutral composition. The thermal expansion in N^, most pronounced 
at high latitudes where the magnetie storm energies are dissipated, causes up- 
welling of air and meridional winds directed toward the equator. The dil'lusion 
velocities induced by these winds upset the static equilibrium for atomic oxygen 
and produce significant variations in its density distribution. 


Our model studies have shown that with an energy source at 100 km associated 
with "hard” auroral precipitation (Rees, private communication), the upward 
transport of O predominates over the loss through meridional winds thus enhanc- 
ing atomic oxygen at high latitudes during magnetic storms -• which is in contrast 
to the observation 

With the energy source at 150 km, presumably due to joule heating (Cole, 1970), 
the meridional wind component is more effective thus reducing atomic oxygen 
at high latitudes and transporting it to low latitudes. This process is most ef- 
fective at altitudes around 200 km and it is consistent with the ionospheric ob- 
servations suggesting a decrease of [Ol / [N^ j during storms (Chandra and 
Hverman, 1969). 

At high altitudes, the thermal expansion in atomic oxygen becomes more 
important and this tends to cancel the wind effects. As a consequence, the 
storm time variations in O are negligible at altitudes above 300 km, which is 
consistent with the observations of Reber, et al. (1970) and Taeusch and Carignan 
(1970). 
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VI. FIGUKK CAPTIONS 


1. I'lic transfer IXmctions, G„, from Volland and Mayr (1971, a) illustrating the 
effieieney for excitation of density components as a function of frequency. 


2. Energy distributions used in computing the wind and temperature fields from 
the model of Volland and Mayr (1971, b). 


3a. Amplitude of the vertical, (V ^ ), and meridional, (V^^ ^ ), wind and temperature, 
(AT), fields computed from the model of Volland and Mayr for the energy 
maximum at 100 km. 


3b. The same parameters as in 3a but for an energy source at 150 km. 

4a. Relative amplitude and phase of the O variations induced by the vertical, 

(V^ 2 )' andmeridional,(V^,^ ), wind fields. The component with subscript, wind, 
shows the combined wind effect. The calculations were made for tm er ergy 
input at 100 km. 

4b. The same parameters as in 4a but for an energy source at 150 km. 

5a. Amplitude and phase of the temperature, (T), and wind, (wind), components and 
the combination of both (without subscript). For comparison, the amplitudes 
and phase of the N 2 concentration are shown. The results correspond to the 
heat input at 100 km. 

5b. The same parameters as in Fig. 5a but for a heat source at 100 km. 

6. The relative variations of O and as functions of latitude at different 

heights. These results agree with observations of the composition and were 
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computed with an energy input that peaks at 150 km suggesting joule dissipation 
as the heating mechanism. 

7. Schematic illustration of the model. 
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Figure 5a. Amplitude and phase of the temperature, (T), and wind, (wind), components and the combination of both (without 
subscript). For comparison, the amplitudes and phase of the N2 concentration are shown. The results correspond to the 
heat input at 100 km. 




temperature, (T)^, and wind, (wind), components and the combination of both (v 
udes and phase of the N2 concentration are shov/n. The results correspond 
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Figure 6. The relative variations of 0 and N 2 os functions of latitude at different heiahts. 
These results agree with observations of the composition and were computed with an 
energy input that peaks at 150 km suggesting Joule heating as the heating mechanism. 
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gure 7. SvSematic illustration of the model. 








